INTEGRALS OF PRODUCTS OF HERMITE FUNCTIONS 
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Abstract. We compute the integrals of products of Hermite functions using the gener- 
ating functions. The precise asymptotics of products of 4 Hermite functions are presented 
below. This estimate is relevant for the corresponding cubic nonlinear equation. 
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Q H . 1. Introduction and statement of the theorem 



In this note, we compute the integrals of products of 4 (normalized) Hermite func- 
tions: 

/oo 
h j (x)h p (x)h q (x)h k (x)dx, j,p,q,k>0 (1) 
-oo 

using the generating functions. The method applies to arbitrary products. The prod- 
uct of 4 Hermite functions is motivated by the cubic nonlinearity in the equation: 

ON 

■9 1 , 3 2 2 N 12 

i— u = -(-^-o + x )u + \u\ u. 



ON 
O 



The special case p = g = was computed in [W], where the author showed stability 
of the harmonic oscillator under the time dependent perturbation: 

V 

V(x, t) = 5\h (x)\ 2 cos(tv k t + <f> k ), (2) 
k=i 

for small 5 and a set of frequencies u> = {ujk} close to full measure. In [W], the precise 
asymptotics: 

W j00k ~ -—=e~W^ forj + /c>l, (3) 
Vj + k 

played an essential role. Using (2) 

ll^/ifclb < -j^ji -> 0, as k -> oo. 

Hence the spatial part of the perturbation diminishes for higher Hermite modes con- 
tributing to stability. 

In this paper, we compute (1) for arbitrary p and q. We prove 
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Theorem. 



r p+i ] 

2 J 1 (A _ U\2 „ 1 ,,_ b ,2 



i _ . V- 1 r(i - *0V 1 (j ~ fc)2 



<?=o 



for all j, k, 



VJTk 

Wj pqk = j + p + q + k odd, 

where [] denotes the integer part and C P;q < a p+q for some a > 1. 

We remark that except for the factor C p ^ q , the estimate in (4) is essentially the same 
as in (3). In particular, the polynomial factor in front of the Gaussian is optimal and 
for fixed p and q, the Schur norm is of the same order as the operator norm. Using 
(4), the result in [W] extends immediately to potentials with exponentially decaying 
Hermite coefficients. The rest of the paper is devoted to the proof of the Theorem. 
We first recall some basic facts about the Hermite functions and the proof of (3) . 

The Hermite functions hj are the eigenfunctions of the harmonic oscillator: 

d 2 



Hh 3 = (-^Z2+ X ) h 3 



de_f ( a t ^ 2 n 

dx" 



with eigenvalues 
and 



A, = 2j + 1, j = 0, 1., 



h j (x) = ^e-^ 1 j = 0,l... (5) 



where Hj(x) is the j th Hermite polynomial, relative to the weight e ^ (H (x) = 1) 
and 

/oo 
e~ x2 HAx)H k (x)dx 
(6) 

So 

1 f°° 

W jpqk = —======== / Hj(x)H p (x)H q (x)H k (x)dx, j,p,q,k>0 (7) 

yj 23 + P + i +k j\p\q\k\ J-OO 

As in [W], the idea is to view the above integral as an L 2 product with the new 
measure e~ 2x and reexpress the products of Hermite polynomials in x as new Hermite 
polynomials in \[2x: 

j+k 

Hj(x)H k (x) =J2 a r H AV2x). (8) 

r=0 
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and similarly 

P+Q 

H p (x)H q (x) = J2beH e (V2x). (9) 

e=o 

Using (8) and (9) in the integral in (7), and assuming (without loss of generality), 
p + q < j + k, we then have 

p+q 



/oo y 1 h 

H j (x)H p (x)H q (x)H k (x)dx = ^aibicg, (10) 
-OO «_n 



1=0 

where 



/oo 
[H £ (V2x)] 2 e- 2x2 dx = 2 i -kl^. (11) 
-oo 

To find the coefficients a r and bg : we use the generating functions as follows. Since 

^ = E h Hn{x) > (12) 

n=0 

oo m 

e 2sx " s2 = £ ^^(x), (13) 
' m! 

m=0 

which can be found in any mathematics handbook (cf. [T] for connections with the 
Mehler formula), multiplying (12, 13), we obtain 

e 2 (*+-)«-(* a +- a ) = ^H n (x)H m (x) 

n\m\ 

n,m 

=e 2(^)V2x-(^) 2 . e -i( t - s)2 (M) 
oo I t+s \£ oo , N 2 r 

= Eff( (^).i^i. E (- ir (^_. 

£=0 r=0 

Using (14) in (8), we note that 

a o = —j+k — —, — 'U + k)\, j + k even 

2 4^(2±^)! (15) 

= otherwise, 

by taking 2r = j + fc, which is the only contributing term. In [W], we computed the 
case p = q = 0: Wjoofc = ao&o c cb which we recall below. 

Lemma. 

"OO 



/oo 
hl(x)hj(x)h k (x)dx 
-oo 



(-l) 1 ^ (j + k)\ [if (16) 

j + k even, 



2i+ k ^]W. (i±*)| V 2 

=0 otherwise. 
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Let 



J = l±±,K j ~ k 



2 ' 

assuming j > k, without loss. When J 1, 



W jk = 



1 + ° { J 



2 ' 



(-l)^J! 



y/2J(J + K)l(J - K)l 



1 Jkl ~VJ 



(17) 

(18) 
(19) 



Proof. (16) follows directly from (15, 7). We only need to obtain the asymptotics in 
(18, 19). This is an exercise in Stirling's formula: 



or its log version 



ji = (iyj2^j(i + — + -J— + ...) 

V V JK 12j 288j 2 ' 



log jl = (j + ^) log j - j + log V2n + ... 



Here it is more convenient to use the latter. Using (17, 20) 

(j + k)l 



log 



2J+fe(i±*)l 



log log 2 2J 



1 



=JlogJ- J + ^logJ + C(J" 1 ). 



So 

using (21). Hence 
WO* 



(j + fc)! 



1 + 



J 



J! 



7rJ 



1 + 



J 



(-l)^J! 



v/2J(J + K)!(J-K)! ! 



J> 1 



which is (18). Using the fact that 



with 



< A, < , for all j > 1, 



(20) 



(21) 



and applying the inequalities (with x = Kj J) : 

<j)(x) d = (1 + x) log(l + x) + (1 - x) log(l - x) > x 2 

for all x G [0, 1) and <j>(x) > ax 2 with a > 1 for x G [7/10, 1), we obtain (19). (When 
x = K/J = 1, (19) follows by a direct computation using Stirling's formula.) □ 
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2. Proof of the theorem 



From (14), the P s term in the RHS is among the terms: 

where I ~ j + k means that I has the same parity as j + k. Equating the coefficients 
in front of the P s term in the second line of (14) and (22), we obtain 

j+k , . j+k-e 

Hj (x)H k (x) =j\k\ H,{V2x). ^ > ° 

e=o,e~j+k z 2 M 2 >■ 

t 



from (8). Similarly 

p+q / i \ P+9- 

H p {x)H q (x) = p\q\ Hi(y/2z) ■ f 



r 
r=0 



i=o,i~ P+ , 2^,!(£±p)! 

E(-i) 3 - m crc^ 

= ^2biHi(y/2x) 



m=0 



from (9). 
So 



( 1) 2 J-^- V~V_-| \k-r r r r k-r 
Z 2 ^ j )■ r=0 

for £ = j + k,j + k- 2, ...(£> 0), 
otherwise, 



and 



for e = p + q,p + q-2,...(£> 0), 

= otherwise. 
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(23) 



(24) 



(25) 



(26) 



Using (25, 26) in (10, 7), assuming p + q + j + k even, otherwise I = 0, we then have 



JUkWah/z p+q ( -\\ j+h t P+q e 

VV jpqk- 2j+p+q + k 2^ ( j + k-£ V( p+q-£ Vp ^2^ [ L) ^e^j + k-i) 



£=0, £~p+q 

I 

( J] (-l) fl - m CTCK^). J+P + Q + k even. 

(27) 



m=0 

For each £ = p + q,p + q — 2, (£ > 0), we then need to estimate 



T (i) _ Vj^kl 1 V~V -\\k-r r'r /~ik—r /no\ 

V 2 /• r=0 

and 

« = #7HE^ E (-1)*-^^- ( 29 ) 

V 2 /• m=0 

Then 

We check that when p + q = 0, the sum in (30) reduces to the term £ = and 
w U + k)\ , . fc 



even, 



same as in (16). 
Estimates on 7jjP. 

We first look at (28). When £ > 0, we need to perform the sum over r. Rewrite 

W = yffl 1 y^, , fc _ r £1 (j + fc-l)! 

jfc 2J+fc (j+fc-i); J r!(£-r)! (fc - r)!(j - + r)! 

_ (~l) fc + I! >p fcl jl 

2j+fc ('+*-*)! ^fW-^ r\(k - r)\ {£ - r)\{j - i + r)\ 

= (-l) k F 1 F 2 , £ = p + q,p + q - 2, ...{£ > 0), 

where F 2 denotes the sum. 
We note that 

(o) i a + k)\ i (j-fc) 2 

7', ' = U , 7 rsj P WI 

Jfc 2^+ fc v /j!H (2±*)| VJT^ 
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(31) 



from the Lemma. So we write 



F = T (o) v_ u + k-e-iy.\ 
' 1 Jk 2I (j + k-i)\\ 



I p I < I r(°) 



jkl 2I O' + *0 [ * ] " 

F 2 can be written as 

k\ j\ 



F 2 = 



r=0 



(32) 



r!(fc-r)! (£ - r)\(j - i + r)\ 

= E(- 1 ) r ^i" r ( 33 ) 

[1/2] 

= £(-l) r CJC£:?\ j>k,j + k>£, 

r=0 

since (t — s) fc (t + s)- 7 = (t 2 — s 2 ) k (t + s) J_fc and both expressions in (33) give the 
coefficients in front of the s e P +h ~ e term (£ < j + k). So 



[e/2] 



w < E 



!-2r 



fc r (j ~ k) 
£ r! (£-2r)! 



We also need to estimate IpJ , I = p + q, p + q — 2, ...{I > 0). Since p + q < j + k, we 
use norm estimates. Comparing (29) with (26) , we have 

where 6^ is as defined in (24). From (24) 

b t J H}{V2x)e- 2x2 dx = J H p (x)H q (x)H £ (V2x)e- 2x2 dx 

< [ J H 2 (x)H 2 (x)e- 2x2 dx} 1/2 [ J H 2 (V2x)e- 2x2 dx] 1/2 



So 



where we used the normalization conditions in (5, 6) and the L°° estimate [T] 

IIMl- < ^712 < a 
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Using (35) in (34), we then have 




So the terms in the sum in (30) can be estimated as follows: 
_ I l-W /-(')! 

-l^fc l e V ^-Z^ ^^JTk> { j + k } r\(£-2r)l 



am 1 . (36) 
<|/ ( ? ) |e^ V - • X £ " 2r - 

< ^__^ , i__ 

where X = . Using Stirling's formula to relate \fl\ and (£ — 2r)! to |! and (| — r)!, 
we have 

(36)<-£=Ve-*-^ ^L, 

1 ^VJ+^o (|-r)! r!(|-r)!' 

for some C > 1. Summing over r and then £, we obtain the Theorem. □ 
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